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ABSTRACT 

We performed a set of cosmological simulations of major mergers in galaxy clusters, in order to study 
the evolution of merger shocks and the subsequent injection of turbulence in the post-shock region 
and in the intra-cluster medium (ICM). The computations have been performed with the grid-based, 
adaptive mesh refinement (AMR) hydrodynamical code Enzo, using a refinement criterion especially 
designed for refining turbulent flows in the vicinity of shocks. When a major merger event occurs, 
a substantial amount of turbulence energy is injected in the ICM of the newly formed cluster. Our 
simulations show that the shock launched after a major merger develops an ellipsoidal shape and gets 
broken by the interaction with the filamentary cosmic web around the merging cluster. The size of 
the post-shock region along the direction of shock propagation is of the order of 300 kpc and the 
turbulent velocity dispersion in this region is larger than 100 km s^^. We performed a scaling analysis 
of the turbulence energy within our cluster sample. The best fit for the scaling of the turbulence energy 
with the cluster mass is consistent with Af^/"^, which is also the scaling law for the thermal energy 
in the self-similar cluster model. This clearly indicates the close relation between virialization and 
injection of turbulence in the cluster evolution. As for the turbulence in the cluster core, we found that 
within 2 Gyr after the major merger (the timescale for the shock propagation in the ICM) , the ratio 
of the turbulent to total pressure is larger than 10%, and after about 4 Gyr it is still larger than 5%, a 
typical value for nearly relaxed clusters. Turbulence at the cluster center is thus sustained for several 
Gigayears, which is substantially longer than typically assumed in the turbulent re-acceleration models, 
invoked for explaining the statistics of observed radio halos. Striking similarities in the morphology 
and other physical parameters between our simulations and the 'symmetrical radio relics' found at 
the periphery of the merging cluster Abell 3376 are finally discussed. In particular, the interaction 
between the merger shock and the filaments surrounding the cluster could explain the presence of 
'notch-like' features at the edges of the double relics. 

Subject headings: hydrodynamics - methods: numerical - galaxies: clusters: general - shock waves - 
turbulence 



1. INTRODUCTION 

Galaxy groups and clusters are the largest virialized 
objects to have arisen in the process of cosmic struc- 
ture formation. These structures have been forming 
by assembling through hierarchical clustering of matter 
driven by gravity. The large scale distribution of mat- 
ter in the Universe has a web-like structure consisting of 
an interconnected network of large filaments, voids, and 
sheets, where cl usters form at the nodes of this mat- 
ter distribution (iBond et all 119961: lEinasto et all 119971: 
iDoroshkevich et aLlUQQlT How such complex structures 
arise out of the primeval density fluctuations has been an 
enduring quest. The cosmic shock waves which are pro- 
duced by accretion and halo mergers play an important 
role in the process of hierarchical structure formation, 
yet their origin and evolution is not very well understood 
(|Bvkov et al.ll2Q08l ). 

Previous numerical simulations fMiniati et al.' '20001; 
iMiniatil l2002t iRvu et a l. 2003; Pfrommer ct al. 2009 
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iSkillman et"all 120081: iVazza et al1l2009[) make a distinc- 
tion between two classes of shocks, namely the 'inter- 
nal' and 'external' shocks. The former group indicates 
weaker (Mach number < 5) shocks found in the 
hot ICM resulting from merging events, whereas the 
latter ones are generated by the in-fall of cold, unpro- 
cessed baryons on accreting structures, thus producing 
shocks with stronger temperature gradients and conse- 
quently larger Mach numbers (A4 > 10). Both exter- 
nal and internal shocks dissipates a large fraction of ki- 
netic energy in the ICM. Strong coUisionless astrophys- 
ical shocks are notably capable of producing high en- 
ergy cosmic-ray par ticles (CR) via diffusive sh ock accel- 
eration mechanism (|Blandford fc Eichleiill987l for a re- 
view). These shocks thus have been proposed as accelera- 
tion sites of the CR particles emittingthe obse rve d non- 
thermal radio emission in clusters jEnsslin ct al. 1998|; 
iRoettiger et al.lll999l: [Miniati et al. |[2001: Mimati .2003i). 
The cosmic 7-ray background may also be generated by 
the same electron populations through inverse-Compton 
scatt ering of the Cosni ic Microwave Backg round Radia- 
tion fLoeb fc Waxman'l 200dlMiniatil[2002tlBagchi et al.l 
,2006; Miniati et af. 2007^ 

Major cluster mergers, where the mass ratio of the in- 
falling halos approaches unity, are among the most en- 
ergetic events of the Universe (jRicker fc Saraziiil 120011 : 
the energy release from a merger of two clusters of mass 
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10^^ Mq, moving with ~ 1000 km s^^ relative to each 
other, may approach 10^* erg). 

Merger shocks are particularly interesting for the 
in.iection of volume-fillin g turbulence in the ICM 
([Subramanian et al.ll2006() . If a turbulent flow is estab- 
lished in the ICM, turbulent dissipation acts on a signif- 
icantly longer timescale than shocks. Turbulence can in 
principle re-accelerate the ambient electrons and amplify 
magnetic fields by the dynamo action. 

In the last few decades, a substantial amount of ob- 
servational evidences for major merger events in clus- 
ters have been collected, mainly in X-rays as well as in 
optical observations, and only recently in radio wave- 
lengths. The radio structures found in galaxy clusters 
are mainly of two types, the central r adio halos and the 
peripheral relics (jFerrari et al.l l2008| ) . Among the pe- 
ripheral radio relics, a special type of 'symmetrical radio 
str ucture' has been found (also known as 'Radio Gischts', 
cf . iKempner et al.l I2004D . These structures are mainly 
found in cluster of galaxies where supposedly a massive 
merger is ongoing. They are thought to be the tracers of 
nearly symmetric shocks propagating through the ICM, 
emanating out of the merging center. Such a structure 
wa s seen for the first time in the cluster A3667, reported 
bv iRottgering et all ()1997D . A spectacular double-radio 
arc which could be part of a gigantic (size of the order of 
2 Mpc) ring-like formation was discovered in the cluster 
Abeh 3376 (A3376) bv Ba gchi et^ (|2006l ). Apart from 
these two, only recently a few more examples of such 
structures have been observed; the filamentary merg- 
ing cluster ZwCl 2341.1-1-000 which shows both diffuse 
radio emission along t he filaments (|Bagchi et al.l 120021 : 
iGiovannini et al.ll2010D an d peripheral double relics sep- 
arated by about 2.2 Mpc (jvan Weeren et all l2009 f) and 
the double r adio relics in th e clust ers A2345 and A1240, 
reported by lBonafede et al.l (|2009f ). 

Concerning the CR acceleration occurring in such 
shocks, two processes are important in the present con- 
text, namely the first-order Fermi acceleration at shock 
front (lDrury iri983: Blandford & Eichlcr 1987) and, per- 
haps, the second-order Fermi acceleration in the turbu- 
lent post-shock region (Schlickeiser .1989) . 

Due to the inherent complexity of the merger events, 
hydrodynamical simulations have been the main inves- 
tigation tool for the study of hierarchical buildup of 
cosmic structures. Furthermore, idealized halo merger 
setups have often been used, in order to simplify the 
physical description of mergers, and to exp lore the pa- 
rameter s pace in a rat h er control l ed way (iHeinz et al. 
2003; As ai et al.l l200l iTakizawal l2005ailbl : lAsai et al. 
200a 120071: iXiang et al.l l2007i r~Clearlv. hydrodynam- 
ical simulations of cluster evolution in a cosmological 
framework are somehow complementary to the above 
cited approach, but provide the final test-bed for the 
idealized studies, and the necessary link to the observa- 
tions. The increase of the availability of computational 
resources and of 'simulation catalogs' make cosmological 
simulations more suitable for a thorough study of merger 
events. 

In this work we have simulated several major merger 
events in a no n-artificial setup wit h the hydrodynami- 
cal code Enzo (|0'Shea et al.|[2005aD . The aim is mainly 
to account for the effect of the propagation of merger 
shocks and their role in injecting turbulence in the ICM, 



in particular in the post-shock region. Emphasis will be 
therefore put on quantitative and morphological features 
of merger shocks and the subsequently generated turbu- 
lence in ICM. In a second part we also analyze the tur- 
bulent flow in the cluster core, a region which is better 
probed by observations, and where the radio halos are 
considered manifestations of the merging activity. Both 
for the turbulence in the post-shock region and in the 
cluster core, the scaling properties with the cluster mass 
have been studied. 

In the present work, we do not explore a detailed con- 
nection between simulations and cluster radio observa- 
tions, although we will be presenting some tentative com- 
parison of our simulated turbulent shock structures with 
the interesting Mpc-scale radio structures observed in 
merging clusters (specifically A3376), revealing appeal- 
ing close similarities. In a future publication we wish to 
present a more detailed comparative analysis of the shock 
induced synchrotron radio structures simulated numeri- 
cally, with the giant peripheral radio structures actually 
observed in dynamically active clusters. 

The work is structured as follows: details of the per- 
formed simulations and the criteria for selecting mergers 
are described in Section[2l The analysis of the performed 
runs and the related results are presented in Section [31 
The robustness of these results with respect to numerical 
convergence is shown in Section |4l and finally the results 
are discussed, compared with observations and summa- 
rized in Section [5l 

2. NUMERICAL SIMULATIONS 
2.1. Setup of the simulations 

We performed a series of cosmological simulations of 
cluster mergers using the Adaptive Mesh Refinement 
(AMR), grid-based h ybrid (N-body plus hydrodynami- 
cal) code Enzo v. 1.0 ()0'Shea et alJl2005aD . Our simula- 
tions assume a flat ACDM background cosmology with 
parameters JIa = 0.7, ilm = 0.3, 57b = 0.04, h = 0.7, 
(Tg = 0-9; and n = I. The simulatio ns have been ini- 
tialized at redshift z = 60 using the lEisenstein fc Hul 
(1999) transfer function, and evolved to z = 0. An ideal 
equation of state was used for the gas, with 7 = 5/3. 
Cooling physics and feedback are neglected, because they 
play only a minor role in the merger problem that we 
investigate in this paper. Furthermore, neglecting cool- 
ing processes imply that the si mulated virialized object s 
obey self-similar scaling laws (i Brvan fc Normanlll998D . 
This turns out useful as it allows us to extrapolate re- 
sults obtained for the modest mass simulated objects 
(M ~ 10^** -^^o) to the more massive and actually ob- 
served structures. 

The simulation box has a comoving size of 
128 Mpc h^^. It is resolved with a root grid (AMR 
level Z = 0) of 64^ cells and 64^ N-body particles. The 
mass of each particle in this grid is 8.3 x 10^"'^ Mq. A 
static grid (/ — 1) is nested inside the root grid. It has 
a size of 64 Mpc h^^ and is resolved in 64^ cells and 
64'^ particles (particle mass 1.03 x 10^^ Mq). Inside this 
grid, in a volume with side of 32 Mpc h~^, a further 
nested grid at Z = 2 is added (64'^ cells and 64'^ particles, 
with mass 1.3 x 10^° Mq), and grid refinement from level 
Z = 3 to / = 6 is enabled. The linear refinement factor 
N is set to 2, allowing an effective spatial resolution of 
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31.25 kpc h-^ at the maximum refinement level. The 
dependence of our problem on resolution issues is dis- 
cussed in Section m The two static grids and the region 
where AMR is allowed are nested on the location of a 
cluster merger, identified according to the prescriptions 
discussed in Section [521 

The mesh refinement is a particularly critical point of 
the setup of such merger simulations. As will be shown 
in the following, a shock is launched soon after a merger 
event, and it propagates from the cluster center to the 
outskirts up to the virial radius and beyond, through a 
medium with gradually decreasing density. An AMR cri- 
terion based only on baryon and dark matter (DM) over- 
density could not track the evolution of this transient in 
a satisfactory way. Therefore, for a better resolution of 
the post-shock regions, we coupled the AMR on over- 
density with a refinement criterion suitable for refining 
the development of turbulence in the vicinity of propa- 
gating shocks. In particular, since shocks are associated 
with flows with negative divergence of the velocity field, 
we used an AMR criterion based on the local variabil- 
ity of the rate of compression of the flow (the negative 
time derivative of the d ivergence d = V • v). This cri- 
terion was developed by ;Sdhmidt_et al. (2009) and used 
in the AMR simulations of llapichino et al. (|2008[ ) and 
llapichino fc Niemeveii (|2008[ ). According to it, a cell is 
considered for refining if the local value of the compres- 
sion rate c(x, t) fulfills the criterion 

ci^,t)>{c),{t)+aX,it) (1) 

where A.; is the maximum between the average (c) and 
the standard deviation of c in the grid patch i, and a is 
a tunable parameter, set to 5.0 for balancing resolution 
and computational efficiency. 
In the over-density criterion, a cell is refined if 

Pi > fiPoniN' , (2) 

where po = 3Hq/8ttG is the critical density. The for- 
mula holds for both baryons and DM. In this work the 
par ameters for over-density are set to /b — /dm = 4.0 
(cf. llapichino fc Niemevedl2008l) . 

2.2. Cluster mergers 

A preliminary, DM-only run was performed in order 
to identify DM halos, and to select merger events to be 
re-simulated with the setup described in the previous sec- 
tion. Th e halos have been ident ified using the HOP al- 
gorithm (jEisenstein fc Hutl 119981 ). 

We focus our study on merger events which occurred 
at z < 0.7 and between halos with mass M > 10^^ Mq 
at the time of merger. As a criterion for selecting ma- 
jor mergers, we rejected events with a mass ratio be- 
tween the merging clumps < 0.5, a value taken as 
the minimum possible ratio at which the cores of the 
two merging halos get destroy ed to form a new cluster 
core ()Salvador-Sole et al]|1998() . Also events starting at 
z < 0.25 were excluded, because the merger shock did 
not have enough time to propagate through the ICM be- 
fore z = in those cases. 

After checking from the mass history that a major 
merger was ongoing, we visually inspected the data for 
prominent merging shocks; in this way, seven represen- 
tative mergers were chosen in our DM-only simulation. 



TABLE 1 

Overview of the performed simulations 



Run 


Keasniit or 
the merger 


[10" Mq] 


M2 
[10" Mq] 


A 


Mf 
[1014 Mq] 


A 


0.7 


2.38 


1.65 


0.69 


1.08 


B 


0.7 


9.75 


5.95 


0.61 


2.21 


C 


0.5 


5.51 


3.29 


0.60 


1.69 


D 


0.5 


2.98 


2.84 


0.95 


0.88 


E 


0.4 


9.86 


5.63 


0.57 


4.73 


F 


0.3 


7.23 


4.23 


0.59 


2.62 


G 


0.25 


2.11 


1.05 


0.50 


0.59 


H 










4.80 



Note. — The masses of the merging clumps (third and fourth 
column) are computed at the merging time (redshift indicated in 
the second column). The fifth column contains the corresponding 
mass ratio for the clumps, whilst in the sixth column the final mass 
of the cluster at 2 = is reported. 



Clusters are then re-simulated as in Section 12.11 Since 
the computational box has periodic boundary conditions, 
the same volume is shifted and the nested grids are cen- 
tered for each run on different locations, corresponding 
to different forming clusters. The selected events span 
different mass ratios, merger redshifts, and total clus- 
ter masses. The features of the simulated mergers are 
summarized in Table [TJ where the runs are identified by 
letters from A to H . Run H is a, special comparison 
case representing a relaxed cluster, which did not experi- 
ence any massive merger in its recent evolution to z = 
(cf. Figure [lOjj). 

3. RESULTS 

3.1. Morphological evolution of the merger shock 

As already stated, our attention will be initially fo- 
cused on the propagation of merger shocks, and on the 
evolution of turbulence in the downstream region. It is 
instructive to introduce those results with a description 
of the main features of the major mergers. Among all 
the simulated merging clusters in Table [TJ the merger F 
has the best morphology with fully developed and nicely 
resolved shock structures, which even propagate beyond 
the virial radius before fading away. This case was there- 
fore chosen as the reference run for our analysis. 

The morphological evolution of merger F is shown in 
Figure [TJ Two sub-clumps approach each other with a 
relative velocity of 980 km s^^, collided for the first time 
at z ~ 0.3 (Figure [T]:) and then pass through a core oscil- 
lation phase. The merging cores are still distinctly visible 
at least until z = 0.1 (Figure [TJ;) before the final coales- 
cence at z = 0.05. The web-like network of filaments is 
also clearly visible around the forming structure. 

The most prominent effect during a cluster merger is 
the evolution of the baryonic component, whose energy 
budget is significantly altered by the event. The gas in 
the ICM is severely attracted in the forming potential 
well, eventually generating a shock wave which propa- 
gates through the intra-cluster gas of the newly formed 
cluster. Part of the kinetic and gravitational energy of 
the merger event is thus dissipated into the ICM. 

This crucial feature of mergers is better followed by the 
evolution of temperature (Figure [5]). The temperature 
increase is first driven by compression at the center of 
the forming cluster (panels [2]2 to [2];), and subsequently 
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z = 0.10 



(g) z = 0.05 



Fig. 1. — The evolution of merger F is shown in density slices. The redshift is indicated at the lower left of each panel, and a identification 
letter is at the lower right. Each panel has a size of 7.7 X 7.7 Mpc and is parallel to the yz plane. The baryon density is color coded, 
and also represented by contours. 

and to the right with respect to the cluster center, by ap- 
plying the Rankine-Hugoniot jump conditions. Between 
z = 0.25 and 0.15, Ai is in the range between 2.5 to 7. 
The larger values are reached for the arc on the right- 
hand side of the panels in Figure [21 resulting in a mild 
asymmetry. 

As this shock propagates out of the ICM of the newly 
formed cluster, it interacts with the surrounding fila- 
ments. The interaction with the web-like cosmic struc- 
ture causes the breaking of the merger shock in separate 
sections, as clearly visible by comparing the bottom row 
of Figure [2] and Figure [TJ This interesting feature is obvi- 
ously not modeled s o far in simulations of idealized merg- 
ers (e.g.lRic kcr & S arazinll200lHRitchie fc Thomasll2m 
[Mitchell et al.ji2009( ). where symmetric bow-like shocks 
propagates unimpeded outwards. 



the shock is launched and propagates outwards. It can 
be followed in the simulation even after it covers few 
virial radii from the cluster center (at z = 0, iivir = 
0.92 Mpc for the cluster in run F). The maximum 
temperature in the central region of the cluster exceeds 
10* K at z = 0.2, several times larger than the cluster 
virial temperature (about 3 x 10'' K). 

The shape of the emerging shock depends on the mass 
of the merging clumps and on the geometry of the 
merger. In the case of merger F, the shock front has a 
roughly ellipsoidal shape, with the arcs more pronounced 
along the merger axis. The morphology of the evolving 
shock front from different lines of sight is shown in Fig- 
ure [21 The propagation velocity of the shock is initially 
up to 1500 km s^^ and only slightly decreases in time 
with the expansion of the shock. We computed the Mach 
number of the two most prominent shock arcs, to the left 
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z=0.25 (d) z=0.20 (e) z=0.15 (f) 




z=0.10 (g) z=0.05 (h) z=0.00 (i) 



Fig. 2. — Same as Figure[T] but temperature is here color coded, with density contours overlayed. In order to display the lower temperature 
structures the upper limit of the temperature (color scale) has been cutoff at 5.0 X lO'^ K (Deep brown-red). The maximum temperature 
at the central region of the cluster goes upto 10* K in certain cases and temperature behind the shock goes beyond 5.0 X 10^ K. 




Fig. 3. — Evolution of the shock in temperature, as seen from slices in three different planes of the computational volume. The panels 
refer to 2 = 0.1 and to slices perpendicular to the x-axis (left column), to the y-axis (central column) and to the z-axis (right column), 
respectively. Each panel has a size of 7.7 X 7.7 Mpc and is cut along the center of mass of the system. 
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3.2. Generation and evolution of turbulence past the 
merger shock 

The development and propagation of a large scale 
shock is a remarkable consequence of the hierarchical 
growth of halos undergoing major mergers. Besides the 
already described effects for the morphology, the ener- 
getics and the heating of the ICM, in this section we 
discuss the injection and evolution of turbulence in the 
post-shock region. 

Velocity fluctuations are a distinctive feature of turbu- 
lent flows. It is therefore straightforward to relate the 
generation of turbulence with the production of vorticity 
occurring at curv ed shock, e xpressed by the curl of the 
Euler's equation (iLandaufc Li fshitz 1959): 



diiJ 

dt ^ ' 



Vp xVp 



(3) 



where the second term at the right-hand side is non- 
vanishing in curved shocks. 

This theoretical expectation is confirmed by our sim- 
ulations. We computed the vorticity from the velocity 
field, using a four-point symmetric method for computing 
the spatial derivatives, with fourth-order space accuracy. 
In the reference run _F, vorticity is produced after the 
merger, just behind the shock, and propagates along with 
it. This process is shown in Figure 21 some level of vor- 
ticity is associated to both clumps before merging (Fig- 
ure |3Ji) and also to the center of the newly formed cluster 
(Figure llf), but the magnitude of u) tracks markedly the 
generation of turbulence in the post-shock region (for 
z < 0.2 , Figures!?^ t olHi). Different from minor mergers 
(ISubran ianian et aT] 120061 : llapichino fc Niemeveii 120081 : 
iMaier et al.ii2OO90 , in the case under consideration a sim- 
ple visual inspection of Figure|3]suggest that major merg- 
ers stir the ICM effectively, resulting in a very volume- 
filling production of turbulence. In Figure HI we also 
observe that values to 5 x 10""'^''' s~"^ corre spond, ac- 
cording to the definitions of IKang et al.l (j2007f) , to a vor- 
ticity parameter (representing the number of local eddy 
turnovers) r > 10, namely to a relatively large vorticity 
and to a fu l l deve loped turbulence. In agreement with 
IKang et al.l ()2007l ). such high values are reached in the 
ICM after the merger, as well as in the post-shock region. 

The size of the post-shock region along the direction of 
propagation of the shock can be roughly estimated from 
Figure [5j In Figure [5] this issue is addressed in a more 
quantitative way, with the analysis of a one-dimensional 
temperature profile computed along a line crossing the 
shock, drawn in Figure [S] Along this line the tempera- 
ture rises abruptly from the pre-shock to the post-shock 
region (from right to left along the x-axis in FigurelHl) and 
then it decreases more gently towards the cluster center. 
As an estimate for the size of the post-shock region, we 
computed its width at half maximum along the line used 
for the profile in Figure [51 given Tmax — 4.3 x 10^ K and 
T,„in ~ 1.3 X 10'' K, the length of the region in Figure HI 
where T > Tmax - AT/2 (with AT = Tmax - Tmin) is 
about 270 kpc h~^. A similar result can be obtained 
from the analysis of the vorticity. This size can be taken 
as a typical order-of-magnitude estimate for regions past 
propagating shocks in major merger events. 

Above in this section we used the vorticity as an in- 
dicative quantity of the turbulent state of the flow. A 



further diagnostic of the magnitude of the turbulent mo- 
tions which is often used is the root mean square (hence- 
forth rms) velocity, defined on a generic domain as 



(4) 



where the sum is performed on the analysis domain, 
and TOj and Vi are the mass contained and the velocity 
magnitude in the cell j, respectively. The quantity {v) 
is the average value of Vi in the analysis volume; from 
an operational point of view, the calculation of Uj-ms in- 
volves the definition of this average reference velocity, in 
order to distinguish between the bulk velocity fiow and 
the velocity fluctuations (cf. iDolag et al.ll2005( ). In our 
case, (v) and Vj^^ns have been computed on a sphere with 
a diameter of 256 kpc h~^, placed inside the post-shock 
region. The size of this analysis region has been chosen to 
be similar to that of the post-shock region, in order to re- 
solve the latter with the largest possible number of cells. 
The shock arc shown in Figure [5] was chosen for the anal- 
ysis, and, to partly avoid spurious fluctuations of Wrms, 
we averaged the results of many (from three to more than 
ten, depending on the evolving size of the region) anal- 
ysis spheres and computed the standard error to show 
its statistical reliability. The method is then applied to 
every output of the simulation, and the results are plot- 
ted in Figure [71 The time evolution of Wmis shows that 
the velocity dispersion is larger when the shock has just 
emerged, and then it decays with a timescale of about 
1.5 Gyr. 

We make use of the velocity dispersion to deflne the 
turbulent pressure support as the ratio of the turbulent 
pressure (Pturb) to the total pressure (Ptot): 



^turb _ 

ftot fcT/(^mp) +i;2„3/3 



(5) 



where k is the Boltzmann constant, fi — 0.6 is the mean 
molecular weight in a.m.u. and nip is the proton mass. 
The turbulent to total pressure ratio (Figure [Sj) peaks at 
relatively large values, close to 10 %, and then it decays 
in a way similar to frms- In both temporal evolutions, 
at about 2.9 Gyr from the shock emergence, one can no- 
tice a mild increase of the turbulence diagnostics, caused 
by the interaction of the merger shock with the reached 
accretion shock. 

The shock propagation does not stir the ICM only in 
the post-shock region. An interesting by-product, high- 
lighted in Figure [9l is the interaction between merger 
shocks and filaments. It is already known that the 
warm-hot baryons flowing along the filaments results in 
a tu rbulent fiow, when it mixes with th e cluster gas 
dNagai fc Kravtsov|[2003l: IMaier et al.ll2Q09l ). Here, as an 
additional effect, we can see that a substantial level of 
vorticity (comparable to that injected in the post-shock 
region) is generated in the regions past the merger shock 
and surrounding filaments. In Figure [9l such 'collars' can 
be clearly observed for the second and third filament 
starting from twelve-o'clock clockwise, and among the 
two filaments on the left-hand side. The turbulence in- 
jection in these zones is probably related to the shearing 
(Kelvin-Helmholtz) instability between the shocked gas 
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Fig. 4. — Same as Figure [T] but slices of the norm of the vorticity = V X v are shown, with density contours overlayed. The vorticity 
has dimension and is reported in s~^ unit in the gray-scale at the upper left corner. 
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Fig. 5. — Temperature slice at z = 0.1 (cf. Figure Y3a\. with a 
size on a side of 3.2 Mpc . Density contours are superimposed. 
The short line crossing the post-shock region is used for computing 
the one-dimensional profile shown in Figure [6] 

moving outwards and the filament gas flowing inwards. 
The level of vorticity is much larger than the turbulence 
associated with the baroclinic generation at the filament 
accretion shock, in regions far from the clusters (well vis- 
ible for the third filament, in the lowest part of Figure |9]). 



Fig. 6. — Profile of temperature along the line shown in the slice 
of Figure [5] The interpolating line is a /3-spline fitting to the data 
points indicated with stars. The horizontal solid line indicates the 
level at T = Tmax - AT/2 (with AT = Tmax - T^in), with the 
lowest temperature past the shock propagation Tmin — 1.3x10^ K, 
and the meaning of the symbols described in the text. 



In Section [5] we further discuss about the link of this fea- 
ture and radio observations of double relics. 
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Fig. 7. — Temporal evolution of frms in the post-shoek region 
in cluster F. The error bars represent the standard errors of the 
computation of Vims- 
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Fig. 8. — As in Figure [T] but the ratio of turbulent to total 
pressure is here reported. 




Fig. 9. — Slice of 12.8 x 12.8 Mpc h^^, showing the squared 
norm of the vorticity uj'^ for the cluster F at z = 0.05. The vorticity 
squared is color coded, whereas the density is represented in shaded 
gray-scale and contours, in order to better highlight the filaments 
surrounding the cluster. 



3.3. Comparison with a relaxed cluster: turbulence 
injection in the ICM 

For a better understanding we compare the properties 
of the merger in run E with those of a relaxed cluster 
{H , in Table [T]). This merger is different from the run F 
used in the previous analysis and it was chosen primarily 
because it has nearly the same final mass as the cluster 
H at z = Q. 

For both clusters we performed an analysis of 
mass-weighted averages of selected quantities (Fig- 
ure nop , computed inside a spherical volume of diameter 
512 kpc h^^ placed at the cluster center. The centers 
are defined as the locations of the peak of DM density, 
as provided by the HOP tool. 

The different accretion histories of the two clusters are 
clearly distinguishable in Figure llOb . where the tempo- 
ral evolution of the cluster mass is reported. A grad- 
ual increase in mass can be noticed for the cluster H, 
whereas cluster E has an abrupt mass gain starting at 
z = 0.5, consistent with a major merger scenario. The 
minor fiuctuation in the evolution of mass are due to the 
uncertainties of the analysis tool. 

The accretion events in the cluster histories are even 
more apparent in the evolution of Wrms fFigurefTUb). This 
quantity has a steep increase in run E during the major 
merger, reaching an average value of about 800 km s^^, 
and then decreases steadily, halving its maximum value 
in about 2 Gyr but still remaining larger than in the 
relaxed cluster, until z — Q. There is some delay between 
the mass increase and the rise of the velocity dispersion, 
probably linked to the emergence of the merger shock. 
In the relaxed cluster, the overall trend of Wrms is more 
regular, although there are two moderate peaks which 
are very likely related with minor mergers. 

The comparison of the temperature evolution (Fig- 
ure II Ob ) is particularly instructive about the role of 
turbulence in dissipating the kinetic energy injected by 
mergers. The two clusters have nearly the same final 
mass and, as expected from simple scaling and virial ar- 
guments, end up at z = with a similar average temper- 
ature T . Whilst for the relaxed cluster T shows almost 
no evolution during the last 6 Gyr, because at z = 0.7 
the structure is already virializing and relaxing, for the 
merging cluster E the temperature increase is slowing 
down only near 2 = 0. The peak of T at z ~ 0.3 
is connected with the propagation of the merger shock 
through the analysis sphere; this event stirs the ICM 
(Figure [TUf)), and the dissipation of the turbulent mo- 
tions is evidently correlated with the increase of internal 
energy. The temperature increase is not limited to the 
peak at z = 0.3 but extends to later times, with a trend 
reminisce nt of the double s hock model for mergers pre- 
sented bv lMcCarthv et alJ (|2007f l. 

Finally in Figure llOl i the turbulent pressure support 
Ptmh/Ptot (equation [5|) is shown. The turbulent contri- 
bution to the total pressure during the major merger is 
significant, being larger than 10% for about 2 Gyr after 
the merger. 

From Figure [TUl i we infer that the level of turbulent 
pressure support (which is equivalent to the correspond- 
ing ratio between energies) typical for minor merger 
events in otherwise relaxed clusters is of the order of 
5% at most, comparable to the value found in the clus- 
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Fig. 10. — Comparison of cluster E (dashed lines) and cluster H (solid lines), through the temporal evolution of selected quantities. 
With the exception of panel a, the computation is based on mass-weighted averages computed on an analysis sphere with the diameter of 
512 kpc h~^, centered at the cluster center. The x axis shows time, indicated by the corresponding redshift, and spans for 6 Gyr. Panel 
a: evolution of the cluster mass. Panel b: rms velocity. Panel c: temperature. Panel d: ratio of turbulent and total pressure. 



ter core by llapichino fc Niemeveii ()2008[ ) . For simplic- 
ity, hereafter we assume that 5% of turbulent pressure 
support is the threshold between relaxed and perturbed 
cluster cores. In case of cluster E, the whole decay time 
to Pturb/-Ptot = 5% is of the order of 4 Gyr. 

When all clusters of our sample are analyzed using this 
criterion, it turns out that the run is a rather typical 
example of turbulent decay after a cluster merger. For 
few clusters, the support level remains even larger than 
10% for most of the evolution after the major merger (see 
Figure [TT]), probably because of later repeated events. 
The shortest decay time, found in one cluster of our sam- 
ple, is 1 Gyr. The relatively large values for the decay 
timescale for turbulent motions fit well with similar esti- 
mates in idealized cluster mergers fe.g.. iRickcr fc Sarazlnl 
l20?)l . 

3.4. Scaling relations for turbulence 

The scaling of ICM features with the cluster mass is a 
useful and informative tool to study the physics of galaxy 
clusters . Self-similar models of clusters (since the seminal 
work of lKaiseJl 19860 rely on the assumption that gravity 
is the main driver of the evolution of such structures, 
whereas departures from simple scaling laws imply some 
more complicate physics at work ( for example, in the 
case of entropy in cluster cores; see lBorgani et al.l 120081 
for a review). 

The scaling of turbulence energy in clusters has 
been investigated w ith a semi-an alytic approach b; 
ICassano fc Bruncttil (|2005D . and by iVazza et al] (|200i 




6 4 

Lookback time [Gyr] 



Fig. 11. — Time evolution of pressure fraction of turbulence in 
all studied clusters of Table [T] The solid black line is used for the 
relaxed cluster H, rest are the merging clusters as indicated in the 
figure legends. 



using SPH simulations. In both cases a scaling of the 
turbulence energy with the cluster mass i^turb oc M", 
with a ~ 5/3, as expected for the thermal energy in the 
virial case. This is co nsistent with the a ssumption, in- 
voked in the model of ICassano fc Brunett i (2005), that 
the turbulence energy injected in the clusters is propor- 
tional to the internal energy of the forming structure, 
which in turn depends on its mass. 

In this section we present a tentative scaling study, 
performed on our simulation set. Obviously, our rela- 
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Fig. 12. — The turbulence energy i?turb (defined in the text) 
is here plotted against the merger mass Afmcrg, for the merging 
clusters of our sample (Table [T]l. The solid line is the best fit to the 
data points, corresponding to a scaling law Sturb M^'^, whereas 
the dashed line (with a slope of 1.66) is the best fit computed by 
excluding the cluster F (the point at logM = 14.0). 

tively small number of runs hardly covering one order of 
magnitude in cluster mass will result in rather large un- 
certainties on the slope of the fit. But nevertheless it is 
interesting to compare the above cited results with those 
from our sample, consisting only of clusters with ongo- 
ing major mergers. Firstly the turbulence injected in the 
cluster core by major mergers is investigated, and then 
we will focus on the post-shock region of merger shocks. 

We define the merger mass Mmcrg as the cluster 
virial mass (including DM) computed with the utility 
enzo_anyl, assuming an over-density S = 200, and the 
turbulence energy i?turb = 0.5 Mgph w^ms' where Wj-ms 
is the velocity dispersion in an analysis sphere with the 
diameter of 512 kpc centered at the cluster center, 
and Afgph is the baryon mass contained in this sphere. 
Both the centers of the analysis sphere and that one for 
the virial mass computation are defined from the peak 
in DM density. The virial radii i?200 in the cluster sam- 
ple range between 0.43 and 0.95 Mpc h^^ right after the 
shock emerged, with a moderate growth at later times, 
so that the analysis sphere has a size that is always well 
contained within the cluster volume. 

In Figure [12] the scaling properties of these quantities 
are studied. The calculation for each cluster has been 
performed at the time of the merger, as reported in Ta- 
ble [TJ The best-fit to the scaling exponent a provides a 
value of 1.6 ± 0.3. The uncertainty is large, but nonethe- 
le ss the result is in agreemen t wit h the previous stud ies 
of iCassano fc Brunettil (pOOl and lVazza et~al] (pOOl l. 

We are not aware of any previous investigation of scal- 
ing laws for the turbulence injected in the post-shock 
region by merger shocks. In Figure fT3l our results on this 
property are presented; the specific turbulence energy 
eturb is defined as 0.5 x v'^^^, where the velocity disper- 
sion is computed in the post-shock region as described 
in Section [221 Here the fitting scaling law has a slope of 
0.74 ±0.25. 

A calculation of the mass in the post-shock region has 
been avoided, because of the related excessive uncertain- 
ties. Assuming that Mpost-shock oc Mmerg, and defining 
the total turbulence energy as i?turb = M^erg eturb > we 
retrieve also in this case a scaling law with a dependence 

5 /3 

compatible with Mmerg, thus indicating also in this case a 



Fig. 13. — Same as Figure [T2l but the effective turbulence energy 
in the post-shock region eturb is on the j/-axis. The solid line is 
the best fit to the scaling law, with exponent a = 0.74, and the 
dashed line (obtained by excluding the outlier point of cluster C) 
has a = 0.88. 

close link between gravitational infall, virialization, and 
injection of turbulence energy. 

It is immediate to see that both the cluster samples in 
Figure [T51 and, to a lesser extent, in Figure [T^j contain an 
outlier. The outlier in Figure [T^l with Afmcrg — 10^"^ M0, 
is the cluster F, while in Figure [T3|(log(Mi„org/-^0) — 
14.4) it is cluster C. We verified by inspecting the data 
that cluster C underwent an off-axis major merger with 
a relatively large impact parameter, leading to a phase of 
core rotation before the final coalescence. It can be spec- 
ulated that, due to this evolution, the energy transferred 
to the merger shock has been smaller th an in the other 
cases , or released on longer timescales (cf. [Mitchell et al.l 
120091) . As for cluster F, a phase of core oscillation has 
been also observed (Section l3.ip . This feature can poten- 
tially affect the analysis also from a numerical viewpoint, 
since it makes difficult for the analysis tool to univocally 
detect the cluster center, introducing inaccuracies in the 
computation of the virial mass and of the turbulence en- 
ergy. 

Removing the outlier from the best-fit calculation pro- 
vide a scaling exponent of 1.66 ± 0.13 (without cluster 
F) in Figure [ni and 0.88 ± 0.07 (without cluster C) in 
Figure [T51 respectively. 

4. RESOLUTION STUDY 

The choice of a refinement criterion and of its relevant 
thresholds and parameters is a very delicate task in AMR 
simulations. Cosmological simulations of galaxy clusters 
are favored by the clumped behavior of these objects: 
criteria based on baryon or dark matter over-density, 
when properly set, are able to catch the relevant struc- 
tures with a good compromise b etween accuracy and 
saving of computation al resources (jO'Shea et al.ll2005bl : 
iHeitmann et al.ll2008[ ). 

Unfortunately the simulations presented in this work 
are extremely challenging from this point of view, be- 
cause we are mainly interested in the evolution of the 
shocked region, whose size at late times is comparable 
with the cluster size, and that for our purpose should be 
carefully refined. Numerical tests lead to the design of 
the setup presented in Section [^T] In this setup the final 
effective resolution is not particularly high, if compared 
with similar grid-based cluster simulations. We verified 
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that the adopted resolution level is adequate for refining 
the post-shock region in major mergers, as the numerical 
tests presented below will show. 

A key feature of grid-based cosmological simulation is 
the force resolution in the DM section of the code. We 
recall that in the second static nested grid at / = 2, 64^ 
N-body particles are allocated in a volume with size of 
32 Mpc h~^, resulting in a particle mass of 1.3 x 10^" Mq 
for our standard setup. In order to study the conver- 
gence properties of our computations, test simulations 
were performed with 32'^ and 128"^ N-body particles, both 
in the root grid and in the two static grids. Both tests 
gave unsatisfactory results: the latter was computation- 
ally not feasible because of constraints on the available 
memory, whilst the former was so poorly resolved that 
the two clumps did not even show a merger shock. For 
this reason, we had to reduce the explored range of the 
resolution study to two further calculations, with root 
grid resolution from / = to / = 2 of 48^ and 96^ grid 
cells and N-body particles. In the following these runs 
will be labeled as i^48, and F96, respectively. Both runs 
allow for six additional AMR levels, so for both the ef- 
fective resolution changes by a factor 1.5 with respect to 
the reference run F. 

One of the most interesting morphological differences 
is a significant time delay in the emergence of the merger 
shock for different resolutions. In particular, the shock 
is launched at z ~ 0.3 for run F96, at z ~ 0.25 for run 
F, and at z ~ 0.05 for run F48. In this last case, the 
evolution is therefore followed only for a very short time 
before z = 0. The reason for this time delay probably lies 
in the transient nature of merger shocks, which appears 
to depend critically on the morphology of the merging 
substructures at different resolutions. In order to permit 
a better morphological comparison, the time delay has 
been compensated in the panels of Figure [TH A qualita- 
tive analysis of the temperature slices shows that there 
is a good morphological match between runs F and F96, 
whereas run FAS cannot be judged on the same grounds, 
because of its excessive delay. 

A more quantitative comparison between the simula- 
tions with different root grid resolution is provided by 
the analysis of the velocity dispersion Umis in the post- 
shock region. We applied the procedure described in 
Section 13.21 for the calculation of Wrms using a number 
of analysis spheres downstream of the merger shock, and 
we repeated for the three runs under examination. Fig- 
ure ll5l shows a moderate convergence of the values of Wmis 
between the runs F96 and F; run i^48 apparently under- 
estimates the velocity dispersion, although an accurate 
comparison is hindered by the too short shock evolution. 
The size of the post-shock region (Figure ^ is also con- 
verged in runs F96 and F, whereas in i^48 the shock has 
not emerged enough at z = to permit the comparison. 

From this resolution study we conclude that runs F96 
and F give results which are consistent each other, 
whereas F48 under resolves the physical processes in the 
ICM which are relevant in this work. Though F96 is 
still marginally manageable within our available compu- 
tational resources, we decided to work on a root grid 
resolution of 64"^, because it turned out to be more suit- 
able for running a whole series of simulations, as those 
performed in this study. 

Besides spatial resolution, obviously the AMR grid res- 



olution plays a relevant role in properly modeling the 
vorticity production. In our numerical scheme, vorticity 
is not explicitly computed by a separate equation, but is 
derived by the velocity field. Is our maximum AMR level 
sufficient for resolving the latter in a sufficient way? 

In order to address this point, we repeated the run F 
with a maximum AMR level 1 — 8, instead of 6; this run 
will be labeled F8. The effective resolution in the run F8 
reaches 7.8 kpc making this simulation equivalent 
to the computations presented in llapichino fc: Niemeved 
(2008) and Maier ct al. (200i). 

After having verified that most of the post-shock re- 
gion benefits of the improved resolution at / = 7 and 8, 
we computed there the evolution of Urms, shown in Fig- 
ure 1151 In this case, the velocity evolution for run F8 
resembles both runs -F96 and F, with a maximum dis- 
crepancy (about 50%) right after the shock is launched, 
where it is most difficult to perform the calculation of 
Urms- No significant time delay has been found in the 
shock evolution, with respect to run F. Interestingly, 
the pressure ratio Pturb/-Ptot in the post-shock region 
(Figure [T6|) is more robust with respect to the spatial 
resolution, for a subtle numerical reason: the increase 
of Wrms (and of Pturb) Icads to an increase also in the 
numerical dissipation and thus to an increase in temper- 
ature and hydrodynamical pressure, roughly balancing 
the turbulent pressure in the pressure ratio (equation [5]). 
An analysis similar to Figure [IT] has been performed us- 
ing the norm of vorticity as diagnostic, with comparable 
results. 

Finally, we observe that the previous resolution study 
can be interpreted as a verification on our simulations, 
namely it shows that our tool accuratel y represents the 
conceptual model (cf. iCalder et al.ll2002[ ). A further, im- 
portant step would be the validation, i.e. demonstrating 
that our setup appropriately describes the nature of tur- 
bulence in the ICM flow. Our study explicitly focuses on 
the generation of turbulence produced by merger shock, 
and to this aim the simulation setup was defined. On 
the other hand, it is known that other stirring mecha- 
nisms can inject turbulence in the ICM. How well can 
they be modeled within our setup, and how much do 
they contribute to the production of turbulence in merg- 
ing clusters? 

We believe that, by definition, the baroclinic vorticity 
production is the main stirrer in the post-shock region, 
and that our setup can follow it in a reliable way, as de- 
scribed above. By extending this analysis to the cluster 
core, other stirring mechanisms have to be considered. 

Tu r bulence production by A GN outflows (jHeinz et al.l 
l2006t iSiiacki fc Springell2006D has not been included in 
the presented setup, but it is not supposed to play a sig- 
nificant role in unrelaxed, non-cool core clusters. Minor 
mergers are also known to effectively inject turbulence 
in the ICM, but only with a s mall volume filling factor 
and in a rather localized way (ISubram anian et al] 120061 : 
llapichino fc Niemeveii 120081: iMaier et al. 200<| These 
mergers inject turbulence via shearing instabilities in the 
wake of the accreted substructure, and this process is 
probably the most prone to be under resolved in hydro- 
dynamical simulations. 

We checked therefore if the higher core resolution in 
run F8 would result in a more accurate modeling of both 
the vorticity and the turbulence generated during minor 
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Fig. 14. — Temperature slices with density contours overlayed, similar to Figure[2] but each having a size of 10.2 X 10.2 Mpc h~^. The 
evolution of runs F9d, F, and F48 is shown in the upper, middle, and lower panels, respectively. The redshifts in the vertical rows do not 
correspond along the different simulations, because of the time delay effect discussed in the text. 
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Fig. 15. — The temporal evolution of the velocity dispersion Urms 
in the pos t sh ock region is computed, using the method described 
in Section |321 for the runs F96, F, FAS and F8 (see legend). The 
different evolutions have been shifted in time, in order to approx- 
imately compensate for the time delay discussed in the text. The 
initial points have a large intrinsic uncertainty, due to the inabil- 
ity of placing an adequate number of examining spheres when the 
post-shock region has a very limited size. 



Fig. 16. — Time evolution of the pressure ratio (equation O in 
the post-shock region, for the runs F (solid line) and F8 (dashed 
line). 

merger events and, in turn, to a larger turbulent pres- 
sure support in the cluster core. We find that, despite 
the higher resolution of run F8 with respect to run F, the 
evolution of both vorticity and the turbulence is qualita- 
tively the same in the two cases. To illustrate the point 
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Fig. 17. — Time evolution of the pressure ratio (equation [5J in 
the cluster core for the simulations F (solid line) and F8 (dashed 
line). 

in Figure [T7] we compare the evolution of Pturb / ^tot for 
both run F8 and run F (see caption for details). We 
only notice quantitative differences which, at most, are 
at the level of 10 to 20%. In order to further inves- 
tigate the resolution issues, in Figure [18] projections of 
the vorticity norm in the cluster core are compared, to- 
gether with contours of DM density (derived from the N- 
body particles with the cloud-in-cell interpolation tech- 
nique). The comparison is performed at z = 0.1, cor- 
responding to one of the largest discrepancies between 
runs F and F8 in Figure [T71 at a lookback time of about 
1.3 Gyr. From Figure [18] we infer that the differences 
are mainly due to small scale structures which are better 
resolved in the F8 run, as shown clearly from the DM 
density contours. According to this result we conclude 
that turbulence injected by the major merger is the lead- 
ing contribution in the core of this cluster. This class of 
mergers acts mainly on length scales which are reliably 
resolved in our reference setup. Other stirring mecha- 
nisms, although possibly underresolved or not resolved 
at all (e. g., turbulence prod uced in the wakes of cluster 
galaxies: iBFegman fc Da"vidlll989l: IKimll2007[) are not ex- 
pected to play a crucial role in turbulence injection in 
cores of merging clusters. 

5. DISCUSSION AND CONCLUSIONS 

In this work the production of turbulence in cluster ma- 
jor mergers has been studied with hydrodynamical AMR 
simulations, with emphasis on turbulence injection in the 
post-shock region and in the cluster core. Representative 
mergers were selected in the simulated cosmological vol- 
ume, and re-simulated in a setup specially designed for 
resolving the post-shock turbulence. 

The morphological analysis was based on a reference 
computation (Section 13. ip with typical features. Right 
after the two sub-clusters fall together, a merger shock 
is launched from the center of the newly formed struc- 
ture, in a roughly symmetric fashion. The shock and the 
associated compression are the main agents of the con- 
version of gravitational to internal energy in the merg- 
ing process. The values of the shock Mach number Ai 
are smaller than the expectations for external accretion 
shocks, but genera lly larger than the typical values for 
internal shocks (cf. iMiniati et al.ll200Cil ). 

The heating is particularly effective past the shock, in 



a region with a size of about 300 kpc In this region, 
as a byproduct of the merger event, some fraction of the 
energy is injected in the ICM as turbulent motions. As 
reported in Figure [S] the energy contribution of turbu- 
lence in this region is at the level of few percent of eint (at 
a length scale of 256 kpc h~^) for a timescale of 1.5 Gyr. 

As for the turbulence at the cluster core, at a length 
scale of 512 kpc the turbulent pressure support (as 
inferred by the velocity dispersion) is about 20% of the 
total pressure. The turbulent to total pressure ratio re- 
mains larger than 10% on a timescale of 2 Gyr, and above 
the threshold of a nearly relaxed cluster (about 5%) for 
4 Gyr, until z = 0. Strictly related to t his point is 
the re latively long timescale for virialization (jPoole et al.l 
120061) . as deduced by the slow temperature increase no- 
ticeable in Figure [TUb . The energy content of the tur- 
bulent motions in the ICM is comparable with previ- 
ous investigations (iDolag et al.l 120051: iVazza et al.l 



previ- 

lHooS 



llapichino fc Niemeveri 120081 : ILau et al.l I2010D and with 
values inferred from o bservations ("Schuecker eralll200l 
IChurazov et~al|[20p. [2008; Werner et al. 200^. 

The analysis presented above shows that a major 
merger event can affect the whole cluster evolution for 
about 4 Gyr or more, a significant fraction of the clus- 
ter history, and much longer than the shock propagation 
timescale. A consequent and interesting question would 
pertain to any observable imprint that such event leaves 
on the cluster structure. A closely related issue concerns 
the bimodality in the correlation between radio and X- 
ra y luminosity of c lusters sho wing a radio hal o , discu ssed 
bv iBrunetti et al.l ([2007) and lBrunetti et al.l (|2009f ). As 
shown in these references, a complete sample of nearby 
ROSAT flux-hmited clusters with Lx (0.1 - 2.4 keV) > 
5 X lO^'' erg s~^ apparently show bimodality, i.e. only a 
small fraction of these clusters have a detectable radio 
halo at 1.4 GHz and their radio and X-ray luminosities 
are correlated, while the rest of the sample have no de- 
tectable radio halo, and only upper limits. In the light 
of the turbulent re-acceleration theory, this bimodality 
was interpreted as originating from the short timescale 
(of the order of 1 Gyr) of the turbulence driving and 
decay. It is currently unclear how this acceleration sce- 
nario and the observed bimodality could be reconciled 
with the theoretical evidences of long turbulence decay 
timescales, presented in our work. Addressing this issue 
requires a careful account of the drop in acceleration ef- 
ficiency following the turbulence decay, which is beyond 
the scope of this paper. 

A more comprehensive comparison with the radio ob- 
servations would require the coupling of the used nu- 
merical scheme with some model of particle accelera- 
tion (M inia ti 2001, 120071 ) . but this is not addressed 
in this study. Moreover, our simulations setup has 
been especially designed for resolving the merger shocks, 
and therefore the spatial resolution in the cluster core 
(roughly corresponding with the effective resolution) is 
relatively low. We will improve upon this numerical 
shortcoming in future works, although the resolution 
study in Section [3] shows convergence of our results. 

The scaling study in Sect ion 13.41 extends the wo rks 
of iCassano fc Brunettil H2005f ) and IVazza eFall (|200^ to 
focus on the special case of major merger: the results, 
though based on a small simulation sample, are consis- 
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Fig. 18. — Projections of the vorticity norm, with DM density contours overlayed, at 2 = 0.1, for the runs F (left) and _F8 (right). 
The projected volume has a size of (640 kpc h^^)"^. The volume is approxi mate ly centered at the location of maximum DM density, and 
contains the sphere for the analysis in the cluster core, as defined in Section 13.41 



tent with the proportionahty between internal and tur- 
bulence energy, since both scale with the cluster mass 
with a dependence close to M^/^ . This makes sense, be- 
cause it is the same scaling law holding for the cluster 
thermal energy in the self-similar model. In other words, 
thermalization and injection of turbulence are two faces 
of a same physical mechanism, the formation of cosmo- 
logical structure through accretion and mergers. 

The results on the shock propagation are compa- 
rable with similar studies performed w i th id e alized 
simulation setups (Ro ettiger eFall 119931 119961 119971: 
[Ricker & Sarazin 2001). The most relevant difference 
with those setups is the interaction between the outgoing 
shock and the filamentary structure which connects the 
newly formed cluster with the surrounding cosmic web. 
The shock propagation is thus hindered in the direction 
of the filaments, resulting in arc-like shapes in tempera- 
ture and vorticity slices (Figures [2] and |4]). 

Here we point out a striking morphological resem- 
blance between the shock wave structures obtained in our 
simulations and those actually observed at radio wave- 
lengths in a few well known merging clusters. In Fig- 
ure[T9]the complex morphology of merging cluster A3376 
(z = 0.046) observ ed with VLA in ra dio and ROSAT in 
X-ray wavelength (jBagchi et al.l [20061 ) is compared with 
the simulated merging cluster F, showing its projected 
temperature and X-ray emission maps. Our aim in this 
work was not of achieving a detailed match. Never- 
theless, a Mpc-scale projected elliptical radio emission 
structure visible in A3376 is remarkably similar to the 
ellipsoidal merger shock front morphology obtained in 
our simulation (Figure [TO]) . Noticeably, the observed ra- 
dio emission and the simulated shock fronts both have 
double arc-like features, that are located at the cluster 
periphery with their concave sides facing the cluster cen- 
ter. Both the observed and the simulated cluster have a 
comparable mass of about 0.5 x 10^^ Mq, major axis of 
elliptical structure of about 2 — 3 Mpc virial radius 
of ^ 1.4 Mpc ICM temperature of about 5 x 10^ K, 
length of the shock structures about 1 Mpc and, 
most noteworthy, both have an e xtension behind th e 
shock front of about 300 kpc h''^ (jBagchi et al.ll2006D . 



The radio morphology matches the morphology of the 
simulated structure in temperature and, more impor- 
tantly, in vorticity, which traces turbulence (see also Fig- 
ures 12 H and p . 

An interesting feature in the eastern radio arc (left- 
hand side in Figure [T51) of A3376 is a 'notch-like' struc- 
ture, where the arc apparently bent inwards towards the 
cluster center. It is natural, in the framework of the 
comparison with the simulations proposed above, to re- 
late this observed morphology with the interaction of the 
merger shock with the filaments of the cosmic web, and 
the subsequent injection of turbulence at the interface 
between the post-shock region and the filament. A sim- 
ilar interpretation, obtained by the analysis of our sim- 
ulation, can be applied also to another double radio arc 
with notch-like featu res, in the merging cluster A3667 
(jRottgering et al.l [1997) . 

Of course, a detailed comparison should make use of 
more sophisticated tools for converting the information 
from the hydrodynamical simulations in synthetic obser- 
vations; this has not been performed here, but is left for 
future work. 

The robustness of the presented results has been tested 
in Section against resolution effects. The implemen- 
tation of additional numerical tools could increase the 
consistency of the turbulence modeling. Although grid 
codes are considered superior to SPH in this prob- 
lem ([Agertz et al.l 120071: iMitcheh et alll2009[ ). for both 
schemes it is not possible to spatially resolve the tur- 
bulent cascade down to the dissipation length scale. A 
way to consistently model the role of turbulence at sub- 
grid scales is provided by Large Eddy Simulations (LES). 
In these tools the unresolved scales are coupled to the 
resolved ones (where the hydrodynamical equations are 
solved) by means of a sub-grid scale model for turbulence. 
Examples of this technique in astrophysical problems are 
found in t he simulations of the explosio i i of type la su- 
pernoyae (iNiemever fc Hillebrandti 119951 : iReinecke et al.l 
I2002t [Schmidt et al.ll2006allbf) as well as in other fields 
(iPope et all 1200 8': 'Scan napieco fc Briiggenl [20081) . Re- 
centlv. IMaieFet al. (2003) coupled LES with AMR, de- 
veloping a tool which is suitable for the study of tur- 
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Fig. 19. — Morphological comparison between the observed cluster A3376 (left) and the simulated cluster F (right). Left panel: composite 
map of radio and X-ray emissions from the cluster A3376. The VLA observed 1.4 GHz total intensity contours (yellow) are: 0.12, 0.24, 
0.48, and 1 mjy/beam. (beam: 20" FWHM). The large ellipse shows an elliptical fit to the peripheral giant radio structure. The central 
color image depicts the thermal bremsstrahlung X-ray emission observed by the ROSAT PSPC instrume nt (~ 12 ks exposure, 0.14-2.0 keV 
band). The red circles mark the positions of the two brightest cluster galaxies. The image is taken from B agchi et al.l |[200f>). Reprinted with 
permission from AAAS. Right panel: simulated projection of cluster F at z = 0.05, with a side of 7.68 Mpc and depth 2.56 Mpc . 
The projected X-ray surface brightness (in the range 0.14 — 2.0 keV) is drawn in colors, and the temperature is overlayed as contours. 



bulence generation in strongly clumped media such as 
galaxy cluster. The application of this novel technique 
to major merger simulations is potentially very interest- 
ing, since the turbulence energy content in these events 
is remarkable and needs to be properly accounted for. 
This problem will be addressed in a forthcoming work. 

The computations described in this work were per- 
formed using the Enzo code, developed by the Labora- 
tory for Computational Astrophysics at the University 
of California in San Diego (http://lca.ucsd.eduj . The 



numerical simulations were carried out on the SGI Altix 
4700 HLRB2 of the Leibniz Computing Center in Garch- 
ing (Germany). S.P. thanks the Deutsche Forschungsge- 
meinschaft (DFG) for providing adequate funding for the 
research and collaboration expenses, and the University 
of Wiirzburg Graduate Schools (UWGS) for providing a 
PhD finishing fellowship (STIBET Abschlufistipendium) . 
Thanks to G. Brunetti for useful discussions, and to 
the anonymous referee for the valuable suggestions and 
constructive criticism which contributed to improve this 
work. 
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